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Abstract 

This paper studies a valuation framework for financial contracts subject to reference and 
counterparty default risks with collateralization requirement. We propose a fixed point approach 
to analyze the mark-to-market contract value with counterparty risk provision, and show that it 
is a unique bounded and continuous fixed point via contraction mapping. This leads us to develop 
an accurate iterative numerical scheme for valuation. Specifically, we solve a sequence of linear 
inhomogeneous PDEs, whose solutions converge to the fixed point price function. We apply our 
methodology to compute the bid and ask prices for both defaultable equity and fixed-income 
derivatives, and illustrate the non-trivial effects of counterparty risk, collateralization ratio and 
liquidation convention on the bid-ask spreads. 
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1 Introduction 


Counterparty risk has played an important role during the 2008 financial crisis. According to the 
Bank for International Settlements (BISQ two-thirds of counterparty risk losses during the crisis 
were from counterparty risk adjustments in MtM valuation whereas the rest were due to actual 
defaults. In order to account for the counterparty risk, recent regulatory changes, such as Basel 
III, incorporate the counterparty risk adjustments in the calculation of capital requirement. On the 
other hand, the use of collateral in the derivative market has increased dramatically. According 
to the survey conducted by the International Swaps and Derivatives Association (ISDA) in 201S^J 
the percentage of all trades subject to collateral agreements in the over-the-counter (OTC) market 
increases from 30% in 2003 to 73.7% in 2013. OTC market participants continue to adapt col¬ 
lateralization and counterparty risk adjustments in their MtM valuation methodologies for various 
contracts, including forwards, total return swaps, interest rate swaps and credit default swaps. 

When an OTC market participant trades a financial claim with a counterparty, the participant 
is exposed not only to the price change and default risk of the underlying asset but also to the 
default risk of the counterparty. To reflect the counterparty default risk in MtM valuation, three 
adjustments are calculated in addition to the counterparty-risk free value of the claim. While credit 
valuation adjustment (CVA) accounts for the possibility of the counterparty’s default, debt valua¬ 
tion adjustment (DVA) is calculated to adjust for the participant’s own default risk. In addition, 
collateral interest payments and the cost of borrowing generate funding valuation adjustment (FVA). 
In the industry, the valuation adjustment incorporating CVA, DVA, FVA and collateralization, is 
called total valuation adjustment (XVA)|^] 

In this paper, we study a valuation framework for financial contracts accounting for XVA. We 
consider two current market conventions for price computation. The main difference in the two 
conventions rises in the assumption of the liquidation value - either counterparty risk-free value 
or MtM value with counterparty risk provision - upon default. Brigo et al. (2012) and Brigo and 


Morini (2011) show that the values under the two conventions have significant differences and large 


impacts on net debtors and creditors. 

With counterparty risk provision, the MtM value is defined implicitly via a risk-neutral expecta¬ 
tion. This gives rise to major challenges in analyzing and computing the contract value. We propose 
a novel fixed point approach to analyze the MtM value. Our methodology involves solving a se¬ 
quence of inhomogeneous linear PDEs, whose classical solutions are shown via contraction mapping 


1 See BIS press release at http://www.bis.org/press/pll0601.pdf 

2 Survey available at http://www2.is da.org/funct ional- areas/research/ surveys/margin-surveys/ 
3 The terminology can be found in Carver (20131 and Capponi (20131, among others. 
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arguments to converge to the unique fixed point price function. This approach also motivates us to 
develop an iterative numerical scheme to compute the values of a variety of financial claims under 
different market conventions. 


In related studies, Fujii and Takahashi (2013) incorporate BCVA and under/over collateraliza¬ 
tion, and calculate the MtM value by simulation. Henry-Labordere (2012) approximates the MtM 
value by numerically solving a related nonlinear PDE through simulation of a marked branching 


diffusion, and provides conditions to avoid a “blow-up” of the simulated solution. Burgard and Kjaer 


(2011) also consider a similar nonlinear PDE in their study of hedging strategies inclusive of funding 
costs, and taking into account closeout payments exchanged at the time when either party of the 
contract defaults. In contrast, our fixed point methodology works directly with the price definition 
in terms of a recursive expectation, rather than heuristically stating and solving a nonlinear PDE. 
Our contraction mapping result allows us to solve a series of linear PDE problems with bounded 
classical solutions, and obtain a unique bounded continuous MtM value. 

Our model also provides insight on the bid-ask prices of various financial contracts. The XVA 
is asymmetric for the buyer and the seller. As such, the incorporation of adjustment to unilateral 
or bilateral counterparty risk leads to a non-zero bid-ask spread. In other words, counterparty risk 
reveals itself as a market friction, resulting in a transaction cost for OTC trades. In addition, we 
examine the impact of various parameters such as default rate, recovery rate, collateralization ratio 
and effective collateral interest rate. We find that a higher counterparty default rate and funding 
cost reduce the MtM value, whereas the market participant’s own default rate and collateralization 
ratio have positive price effects. For claims with a positive payoff, such as calls and puts, we establish 
a number of price dominance relationships. In particular, when collateral rates are low, the bid-ask 
prices are dominated by the counterparty risk-free value. Moreover, the bid-ask prices decrease 
when we use the MtM value rather than counterparty risk-free value for the liquidation value upon 
default. 

The recent regulatory changes and post-crisis perception of counterparty risk have motivated 


research on the analysis of XVA. Brigo et al. (2014) consider an arbitrage-free valuation framework 


with bilateral counterparty risk and collateral with possible rehypothecation. Capponi (2013) studies 
the arbitrage-free valuation of counterparty risk, and analyzes the impact of default correlation, 


collateral migration frequency and collateral re-hypothecation on the collateralized CYA. Thompson 


(2010) analyzes the effect of counterparty risk on insurance contracts and examines the moral hazard 


of the insurers. Brigo and Chourdakis (2009) focus on the valuation of CDS with counterparty risk 


that is correlated with the reference default. Hull and White (2012) investigate the wrong way risk 
- the additional risk generated by the correlation between the portfolio return and counterparty 
default risk. 

Let us give an outline of the rest of the paper. In Sect. [2j we formulate the MtM valuation of 
a generic financial claim with default risk and counterparty default risks under collateralization. In 
Sect. [3j we provide a fixed point theorem and a recursive algorithm for valuation. In Sect. [4j we 
compute the MtM values of various defaultable equity claims and derive their bid-ask prices. In 
Sect. [5j we apply our model to price a number of defaultable fixed-income claims. Sect. [6] concludes 
the paper, and the Appendix contains a number of longer proofs. 
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2 Model Formulation 


In the background, we fix a probability space (0, J-, Q), where Q is the risk-neutral pricing measure. 
In our model, there are three defaultable parties: a reference entity, a market participant, and a 
counterparty dealer. We denote them respectively as parties 0, 1, and 2. The default time r* of party 
i E {0,1,2} is modeled by the first jump time of an exogenous doubly stochastic Poisson process. 
Precisely, we define 


Tj = inf [t > 0 : 



( 2 . 1 ) 


where {£}}*= 0 . 1,2 are unit exponential random variables that are independent of the intensity pro¬ 
cesses (A^)t>o, i E {0,1, 2}. Throughout, each intensity process is assumed to be of Markovian form 
A« = St, X t ) for some bounded positive function A ^(t,s,x), and is driven by the pre-default 

stock price S and the stochastic factor X satisfying the SDEs 


dS t = ( r(t , X t ) + A(°)(t, S t , X t )) S t dt + aft, S t ) S t dW t , (2.2) 

dX t = b(t,X t )dt + Ti(t,X t )dWt. (2.3) 


Here, (Wt)t>o and (Wt)t >0 are standard Brownian motions under Q with an instantaneous correla¬ 
tion parameter p E (—1,1). The risk-free interest rate is denoted by rt = r(t,Xt ) for some bounded 
positive function. At the default time tq, the stock price will jump to value zero and remain worth¬ 


less afterwards. This “jump-to-default model” for S' is a variation of those by Merton (1976), Carr 


and Linetsky (2006), and Mendoza-Arriaga and Linetsky (2011). 


2.1 Mark-to-Market Value with Counterparty Risk Provision 


A defaultable claim is described by the triplet ( g,h,l ), where g(Sx,XT) is the payoff at maturity 
T, (h(St, Xt))o<t<T is the dividend process, and 1{tq,X T o ) is the payoff at the default time tq of the 
reference entity. We assume continuous collateralization which is a reasonable proxy for the current 


market where daily or intraday margin calls are common (see Fujii and Takahashi 2013). For party 


i E {1, 2}, we denote by <5,; the collateral coverage ratio of the claim’s MtM value. We use the range 
0 < 5i < 120% since dealers usually require over-collateralization up to 120% for credit or equity 


linked notes (see Ramaswamy 2011, Table 1). 


We first consider pricing of a defaultable claim without bilateral counterparty risk. We call this 
value counterparty-risk free (CRF) value. Precisely, the ex-dividend pre-default time t CRF value 
of the defaultable claim with ( g , h, l) is given by 


n(f, s,x) := 


e - ff( r u+\u ) ) du g (s Tl X T ) + f T e -^+^ du (h(S u ,X u ) + X^l{u,X u )^j du 


(2.4) 


The shorthand notation Et sx [-] := E[- \St = s, Xt = x] denotes the conditional expectation under 
Q given St = s, X t = x. 
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Incorporating counterparty risk, we let r = minjro, Ti,T 2 }, which is the first default time among 
the three parties with the intensity function A (t, s, x ) = Ylk=o A®(f, s, x). The corresponding three 
default events (t = To}, {r = t{\ and {r = T 2 } are mutually exclusive. When the reference entity 
defaults ahead of parties 1 and 2, i.e. r = to, the contract is terminated and party 1 receives 
1(to,X tq ) from party 2 at time to- When either the market participant or the counterparty defaults 
first, i.e. t < to, the amount that the remaining party gets depends on unwinding mechanism at 
the default time. We adopt the market convention where the MtM value with counterparty risk 


provision, denoted by P, is used to compute the value upon the participant’s defaults (see Fujii and 


Takahashi 

2013 

Henry-Labordere 

2012) 


Throughout, we use the notations x + = xl{ 3 .> 0 } and x~ = -sl; x< o}. Suppose that party 2 
defaults first, i.e. t = T2. If the MtM value at default is positive ( P T2 > 0), then party 1 incurs a 
loss only if the contract is under-collateralized by party 2 (82 < 1) since the amount 82 P + 2 is secured 
as a collateral. As a result, with the loss rate L 2 (i.e. 1 - recovery rate) for party 2 , the total loss 
of party 1 at T 2 is L 2 (1 — 82 ) + P+ 2 - On the other hand, suppose that the MtM value is negative 
(. P T2 < 0). Party 1 has a loss only if party 1 puts collateral more than the MtM value P T2 , i.e. the 
contract is over-collateralized (eh > 1). In this case, party l’s total loss is the product of the party 
2’s loss rate and the exposure, i.e. L 2 (<5i — 1 ) + P t T- Therefore, the remaining value of the party l’s 
position at the default time T 2 is 


P T2 ~ L 2 (1 - 5 2 ) + P+ - L 2 (5i - 1) + P~ . (2.5) 

Next, we consider the case when party 1 defaults first, i.e. t = ti. We denote by L\ the loss 
rate of party 1. If the MtM value of party l’s position at the default is negative (P Tl < 0) and the 
contract is under-collateralized (5i < 1), party 2’s loss is L\ (1 — 5i) + P~. Similarly, when the MtM 
value is positive (. P Tl > 0) and the contract is over-collateralized (<5i > 1), party 2 incurs a loss of 
the amount L\ (82 — l) -1- ! 3 ^. Because of the bilateral nature of the contract, party 2’s loss is party 
l’s gain. Therefore, at the default time ti, the value of party l’s position is 

P n +Li( 1 - *i) + P~ + U (82 - 1) + P+ • (2.6) 

Moreover, the market participant is exposed to funding cost associated with collateralization 
over the period since the collateral rate and funding rate do not coincide with the risk-free rate. 
When the liquidation value of the contract Pf is positive to party 1 at time t, party 2 posts collateral 
82 T) + to party 1. To keep the collateral, party 1 continuously pays collateral interest at rate C 2 to 
party 2 until any default time or expiry. On the other hand, when Pf is negative to party 1, party 1 
borrows <5i T) - to post collateral to party 2. As a result, party 1 receives interest payments at rate 
c\ proportional to collateral amount. We call c,; the effective collateral rate of party i (i = 1,2), 
which is the nominal collateral rate minus the funding cost rate of party i. The rates ci and C 2 can 
be both negative in practice if the funding costs are high. Therefore, party 1 has the following cash 
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flow generated by the collateral and effective collateral rates: 


!{t<r} (ci <5i P t ~ c 2 5 2 P t + ) , 0 < t < T . 


(2.7) 


The aforementioned cash flow analysis implies that the pre-default MtM value with counterparty 
risk (CR) provision is given by 


P(t,s,x) = E t , 


z -f t T Pu+Xu)du 9 ^ St ^ Xt) + r e -f t u {r v +X v )dv (h(S u ,X u ) + Af *(*,*«)) 


du 


+ j e -i;p v+ x v ) d v ((1 _ L2 (1 _ s 2)+) P + _ (i + l 2 («?! -1) + ) p-) du 

+ [ T A« e-K^+x^dv ((1 + Li (S 2 - 1)+) P+ - (1 - L, (1 - «!)+) P~ ) du 


+ [ e -f t u P,+x v )dv (ci8\P~ - c 2 5 2 P+) du 


( 2 . 8 ) 


The first line accounts for the terminal cash flow, the dividend, and the payoff at the reference 
asset’s default (r = tq). The second line and the third line are the cash flows at party 2’s default 


(r = t 2 ) in ( |2.5[ ) and party l’s default (r = ti) in (2.6), respectively. The last line results from the 
collateral and effective collateral rates in (|2.7|). To simplify, we introduce the following notations 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


f{t , s, x ) = r(t, x ) + A (t, s, x) , 

a(t, s, x) = L 2 A (2) (t, s, x) (1 — <5 2 ) + - Li A^ ( t , s, x) ( 5 2 - 1) + + c 2 8 2 , 

/ B(t , s,x) = Li A (1) (t, s,x) (1 - 5i) + - L -2 A (2) (f, s,x) (<5i - 1) + + ci , 

f(t, s,x,y ) = h(s,x ) + A (0 )(f, s,x) l(t,x ) + (A (1 ^ + A (2) - /3)(t, s, x)y + (/3 - a)(t, s,x) y + . 


This allows to express (2.8) in the equivalent but simplified form: 


P(t, s , cc) — 


e- It du g (g Ti Xt) + e- s: fv S U 1 X U , P u ) du 


(2.13) 


where ft = f(t, St,X t ) as defined in (2.9). 

Remark 2.1. As an alternative of MtM value with CR provision, the liquidation value at the time 
of default can be evaluated as the CRF value of the claim. In other words, at the default time 


r < tq, the liquidation value is evaluated as IR rather than P T . Replacing P u in (2.13) with IR for 


t < u <T gives the MtM value without CR provision (see Henry-Labordere (2012)): 


P{t , S, x) — Ej jSi a 


e" £ du g(S T , X T ) + e~ % dv f(u , S U ,X U , IR) du 


(2.14) 


To conclude this section, we summarize the symbols and their financial meanings in Table [l 
which we will use frequently throughout this paper. 
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Symbol 

Definition 

Symbol 

Definition for party i £ {1,2} 

P 

MtM value with CR provision 

Ri 

Recovery rate 

P 

MtM value without CR provision 

Ci 

Effective collateral rate 

n 

CRF value 

Si 

Collateralization ratio 

TO 

Default time of reference asset 

Ti 

Default time 

A(°) 

Default intensity of reference asset 

A « 

Default intensity 


Table 1: Summary of notations. 


2.2 Bid-Ask Prices 

In OTC trading, market participants, such as dealers, may take a long position as a buyer or a short 
position as a seller. Without counterparty risk, the buyer’s CRF bid price II b {t, s , x) for a claim with 

by replacing 


payoff (g,h,l) is given by (2.4). The MtM value of the seller’s position satisfies 
( g , h, l ) with (—g, —h, —l), the negative of which gives the seller’s CRF ask price II s (t, s, x). In fact, 
the bid-ask prices are identical, i.e. II 6 (f, s, x) = n s (i, s, x). 

Similarly for the case with counterparty risk provision, the buyer’s bid price is P b (t,s,x ) = 


P(t,s,x ) as in (2.13). The seller’s ask price is given by 


P s (t,s,x) = E tiSia 


e“ $ du g(S T , X T ) + / e- -A" */(«, S u , X u , P°) du 


(2.15) 


where 


f(t,s,x,y ) = h(s,x) + A ( °)(t, s,x)l(t, x) + (A (1 ) + A (2 ) - /3)(t,s,x)y — (/3 — a)(t,s,x)y . (2.16) 


Since f(t, s, x , y) is different from f(t , s, x, y ) in (2.12), the symmetry observed in the CRF prices 
generally no longer holds in the presence of bilateral counterparty risk. Most importantly, such an 
asymmetry generates bid-ask spreads for defaultable claims. For any contract with counterparty 
risk provision, the participant can quote two prices: P b (t , s,x) as a buyer or P s (t, s,x) as a seller. 
In addition, since the payoff components ( g , h, l ) can be negative, the bid and/or ask prices also can 
be negative (see Figure [4]). 

The total valuation adjustment (XVA) is defined as a deviation of the MtM value from the CRF 
value, namely, II — P b for a long position and P s — II for a short position. The bid-ask spread 
accounting for the XVA with CR provision is defined as S(t, s, x) = P s (t, s, x) — P b (t , s, x). 


The two factors a and (3 in (2.10) and (2.11) that appear in / and / summarize the effects of 


counterparty risk and collateralization on the bid-ask prices. Specifically, a explains the effect of 
positive counterparty exposure of the MtM value P+ while /3 explains the effect of negative exposure 
P~. When the two parameters have the same value (a = f3), the two functions / and / in (2.12) 


and (2.16) are identical. Therefore, the bid-ask prices P° and P s are equal. Such a price symmetry 


also arises in a number of other scenarios: (i) when both parties have perfect collateralization ratio 
(<5i = §2 = 1) and the same effective collateral rate (ci = C 2 ); (ii) when both parties have zero 
collateralization ratio (<5i = 82 = 0) with the same effective default rate (L\ A^ 1 ) = L 2 A^ 2 -*), and (iii) 
when both parties have the same effective collateral rate (ci = C 2 ) with the same effective default 
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rate and collateralization ratio (L\ A^ 1 * = L 2 (5 1 = 82 )- 

Remark 2.2. When the counterparty risk-free value II is used to estimate the liquidation value upon 
default, the seller’s bid price is given by 


P s (t, s,x ) = E t , s ,x 


e“ ^ du g(Sr, X T ) + / e" K dv f(u, S Ul X u , U u ) du 


(2.17) 


where / is defined in (2.16). In contrast to (2.13), the price function on the LHS does not appear 
on the RHS. 


3 Fixed Point Method 


The defining equation (2.13) has a recursive form whereby the price function P appears on both 


sides. Denote the spacial domain by V := M + x M. For any function w G C&([0, T] x T>, M), we define 
the operator M. by 


(. Mw)(t,s,x ) = E ti , 


e“ & r '“ du g{S T , X T ) + e~ K dv f{u, S U ,X U: w{u , S u , X u )) du 


(3.1) 


Then, we recognize from (2.13) that the MtM value with counterparty risk provision satisfies P = 


MP. This motivates us to show that the operator M. has a unique fixed point, and therefore, 
guarantees the existence and uniqueness of the MtM value P. 


We discuss our fixed point approach by first showing that the operator A4 defined in (3.1) 


preserves boundedness and continuity. To this end, we outline a number of conditions according to 


Heath and Schweizer (2000). 


(Cl) We define 

r(t,s,x) = 


r(t, x) -I- A (°)(t, s, x) 
b(t, x) 


and Tj(t,s,x) = 


a(t,s)s 0 

pri(t,x) y/l- p 2 r](t,x) 


The coefficients T and £ are locally Lipschitz-continuous in s and x, uniformly in t. That is, 
for each compact subset F of T>, there is a constant Kp < 00 such that for -i/> G {T, £}, 

\ip(t,si,xi) - ip(t,s 2 ,x 2 )| < K f \\(si,xi) - (S 2 , £2)11 Vt G [0, T], (si,aq), (s 2 ,x 2 ) G F, 

where || ■ j| is the Euclidean norm in M 2 . 

(C2) For all ( t,s,x ) G [0, T) x V, the solution (S', X) neither explodes nor leaves V before T, i.e. 


sup ||(S n , X u )|| < 00 = 1 and Q( (S u , X u ) G V , Vu G [t, T] ) = 1 . 

, t<u<T / 


(C3) The functions h and g are bounded and continuous, and r, l and A^, i G {0,1, 2}, are positive, 
continuous and bounded. 
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Lemma 3.1. Given any function w £ C&([0, T] x D, M), it follows that v := M.w £ Cf,([0, T] xP,M). 

Proof. The boundedness of v follows directly from that of w , h, g, r, and (see condition (C 3)). 
To prove the continuity of v, we first observe that 

(t, s, x) ^ e -ft T ^ du g(S T , X T ) + ^ e~^^ dv f(u, S u , X u , w(u, S u , X u )) du (3.2) 


is continuous Q-a.s. Indeed, the continuity of (S', A") implies that the mapping (t,s,x) *->• g(Sx,XT) 
is continuous Q-a.s. Also, (t,s,x,u) i->- f(u,S u ,X u ) and ( t,s,x,u ) H > f(u,S u ,X u ,w(u,S u ,X u )) are 
uniformly continuous and bounded Q-a.s. on compact subsets of [0, T]xDx[t, T], Hence, the mapping 
in (3.2) is continuous Q-a.s. Recall that v = Xiiv is the expectation of the RHS in (3.2), which 
is bounded continuous, so v is also continuous for (t, s, x) £ [0, T] x V by Dominated Convergence 
Theorem. □ 


3.1 Contraction Mapping 

Next, we show that the mapping At is a contraction. By the boundedness of a(t, s, x), (3(t , s, x) and 
A «(t, s , x) for i £ { 0 , 1 , 2 }, we can define a finite positive constant by 


L = sup 

(t,s,ir)E[0,T] xT> 


,s,x) + A ( 2 )(t, s,x) - /3(t, s,x)| + \/3(t, s,x) - aft , s,x)||. 


Proposition 3.2. The mapping At defined in (3.1) is a contraction on the space Cf,([0,T] x T>,} 
with respect to the norm 


= 1<yT ^\w(t, s, x)| , 

for L < 7 < oo. In particular, At has a unique fixed point w* £ C^QO, T] x DA 


|u >|| 7 := sup e 

(t,s,i)6[0,T]xP 


(3.3) 


Proof. From (2.12), we observe that \f(t, s, x,y±) — f(t, s, x, ij 2 )| < L \yi — 7 / 2 1, for ( t, s, x) £ [0, T]xD. 
This implies / is Lipschitz-continuous in y, uniformly over (t, s, x). By Lemma 3.1 the operator At 
maps Cb([0,T] x D,M) into itself. For ( t,s,x ) £ [0, T] x V, w\,W 2 £ Cfc([0,T] x D,M), and 7 > 0, 
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we have 


»~7 (T—t) | 


e ‘'|(.A4u;i)(t, s, ®) — (M.W 2 )(t, s, x)| 


= e" 7 ^ 


E 


't,S,X 


e ^ rvdv (f(u,S u ,X u ,wi(u,S u ,X u )-f(u,S u ,X u ,W 2 (u,S u ,X u )))du 


< e _ 7 (T_t) E- 


t,S,X 


(**) frr, ,S 

< e _ 7 (T_t) E- 


t,S,X 


f ( u , *5 W , f {u, S u , X u , w>2 (^5 -X^)) 


e - 7 (r-«) L | ti;i ( u>5ttjjYu) _ W2 (^s n) X u )|e 7 ( T - u )du 


du 


it 

(in') r 

< e _7 ^ T_ ^L||u;i — it72||-y / e 7 ^ T_u ^dit 

L 

< — 101 - 102 U • 

7 


We have used the condition that f„ > 0 for (i), and the fact that / is Lipschitz in y for (n). 


The inequality (m) is implied by the norm in (3.3). As a result, for any 7 > L > 0, M is a 
contraction. S 


The norm 


is equivalent to the supremum norm 


on the space C&([0, T] x T>,M). A 


similar norm is used in Becherer and Schweizer (2005) and Leung and Sircar (2009) in their studies 


of reaction diffusion PDEs arising from indifference pricing. 

Using the fact that Xi is a contraction proved in Proposition 3.2 there exists a sequence of 
functions (p ( n )) n>0 that satisfy p( n+1 ) = XiP( n \ Vn > 0 , and the sequence converges to the fixed 
point P. The convergence does not rely on the choice of the initial function. Indeed, one can simply 
pick any bounded continuous function as a starting point, e.g. P = 0 V(t,s,x), and iterate to 
have a sequence (P( n ')n>o that resides in Cb([0,T] x V, M). 

Furthermore, we can show that for each n > 1, P^ = P^ n \t,s,x) is a classical solution of the 
following inhomogeneous PDE problem: 


dp ( n ) 

dt 


+ CP^ n) - fit , s, x ) P (n) + fit, s, x, P^" 1 )) = 0 , 

P (n) (T,s,x) =g(s,x), 


(3.4) 


where the operator L is defined by 


„ 1 9 2 1 . x2 d 2 , . , . d 2 

c := 7 (t - s) ' s a? + 2 n{t ' x) a^ + n{t ' x)a{t ' a)a aoi 

x 5 ,, . <9 

+ r(t,s,s)s — + 6 (t,x)— . 


(3.5) 


In order to prove the result, we need the following additional conditions, adapted in our notation 


from (A3') — (A3d') of Heath and Schweizer (2000). 


(C4) There exists a sequence (P n ) ne ^ of bounded domains with closure D n C D such that L)^ =1 D n = 
D and each D n has a C 2 -boundary. 
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As in 


Heath and Schweizer 


(2000), one can take V n = [—,n] x [—n,n\ C M+ xl. For each n, we 


require that 

(C5) b(t,x), a(t,s,x ) := X(t, s, x) X*(t, s, x), and r(t,s,x ) be uniformly Lipschitz-continuous on 
[0,T] x D n , where X* denotes the transpose matrix of X, 

(C 6 ) a(t,s,x) be uniformly elliptic on M 2 for ( t,s,x ) E [0, T) x D n , i.e. there is 5 n > 0 such that 
y t a(t,s,x)y > 5 n \\y\\ 2 for all y E M 2 , 

(C7) f(t,s,x,y) be uniformly Holder-continuous on [0,T] x D n x M. 

The conditions (Cl) - (C7) are quite general, and they allow for various models, including the 
Heston, CEV, and thus, geometric Brownian motion models for equity, and the Ornstein-Uhlenbeck 


and Cox-Ingersoll-Ross models for the stochastic factor X (Heath and Schweizer, 2000, Sect. 2). 


The triplet ( g,h,l ), default intensities A® and interest rate r can be easily chosen to satisfy the 
boundedness and continuity conditions in (C3), as we will do in our examples in Sections [4] and [5j 

Theorem 3.3. Under conditions (Cl) — (C7) , there exists a sequence of bounded classical solutions 
(p( n )) c£’ 2 ([0, T)xV, M) of the PD E problem (3.4) that converges to the fixed point P E Cfc([0, T) x 

T>, M) of the operator A4. 


We provide the proof in Appendix | A. 1| The insight of Proposition 3.2 and Theorem 3.3 is that 
we can construct and solve a series of inhomogeneous but linear PDEs whose classical solutions 


converge to a unique fixed point price function P as in (2.13). Recent studies by Burgard and 


Kjaer (2011) and Henry-Labordere (2012) evaluate the MtM value P(t,s,x ) by working with the 
associated nonlinear PDE of the form: 


dP 

— + CP - r(t, s, x) P + f(t, s, x, P) = 0 , 


(3.6) 


for (t, s, x) E [0, T) xT>, with terminal condition P{T , s, x) = g(s, x), for (s, x) E V. The nonlinearity 


of (3.6) poses major challenges on analyzing and numerically solving for P. Henry-Labordere (2012) 


provides a method to approximate the solution that involves replacing the nonlinear term / with a 
polynomial and simulating a marked branching diffusion. This method, however, does not guarantee 
that the solution from simulation will resemble the solution of the nonlinear PDE, and does not 


ensure any regularity, such as continuity or boundedness of, either solution. Henry-Labordere (2012) 


provides conditions on the chosen polynomial to avoid a “blow-up” of the simulation algorithm. 
In contrast, our fixed point methodology circumvents this issue by establishing that the pricing 


definition in (2.13) is a contraction mapping, as opposed to working with the nonlinear PDE. As 


a result, we solve a series of linear PDE problems with bounded classical solutions. In the limit, a 
unique bounded continuous MtM value P is obtained. 


3.2 Numerical Implementation 

Our contraction mapping methodology lends itself to a recursive numerical algorithm. As mentioned 


in the previous section, we iteratively solve a sequence of linear inhomogeneous PDEs (3.4). At 
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each iteration, the error is measured in terms of the maximum difference between two consecutive 
solutions P n ' ) and p( n ~ 1 '> over the entire domain [0, T] x T>. We continue the iteration procedure 
until the error is less than the pre-defined tolerance level e. 

For implementation, we use the standard Crank-Nicolson finite difference method (FDM) to 


obtain the values (see, among others, Wilmott et al. (1995) and Strikwerda (2007)). We restrict 
the domain [0, T\ x T> to a finite domain T> = {(t,s,x) : 0 < t < T, X < x < X,0 < s < 
S}. The parameters S, X_ and X are sufficiently large enough to preserve the accuracy of the 
numerical solutions. We discretize the function p( n \t,s,x ) as P^ n \ti, Sj, Xk) where i G {0, ...,1V}, 
j € {0,...,M} and k G {0, ...,L} with At = T/N, As = S/M , Ax = (X — X)/L and ti = iAt, 
Sj = j As, Xk = k Ax. Our numerical procedure is summarized in Algorithm [lj 

Algorithm 1 Fixed Point Algorithm for Evaluating the MtM Value P 

set n = 1, P + = P 
solve for P^ from PDE (3.4) 
set e = H-PC 1 ) — -P(°)||oo 
while e > e do 

set n = n + 1, P + = pO- 1 ) 
solve for P^ from PDE (3.4) 
set e = ||P( n ) -PMIU 
end while 
return P 


For the CRF value, we solve the linear PDE associated with II = II(f, s,x) in (2.4), namely, 


<9n 

at 


+ c n 


r(t, s, x) n(t, s, x) + h(s, x) + s, x) l(t, x) = 0 , 


(3.7) 


for (t, s , x) G [0, T) x V, with terminal condition II(T, s, x) = g(s, x), for (s, x) G V. The CRF value 
becomes an input to the PDE problem for the MtM value without provision, given by 


d P 
~dt 


+ CP 


fit, s,x) P + f(t, s, x, II(t, s, x)) = 0 , 


(3.8) 


for ( t , s, x) G [0, T ) x V, and P(T, s, x) = g(s, x), for (s, x) G V. Again, we apply the Crank-Nicolson 
FDM method to compute their values. 


Remark 3.4. The assumption on the boundedness of the default intensities does not encapsulate the 
local intensity model, where the reference default intensity function is of the form: A (t, s ) = cs~ p , 


for some p, c > 0. See Carr and Linetsky (2006); Carr and Madan (2010); Linetsky (2006); Madan 


and Unal (1998), among others. As a close alternative, one can cap the exploding intensity and 


set A (t, s ) = cs~ p A M for some arbitrarily chosen large constant M. With this modification, our 
model and contract mapping results still apply. Since the default intensity function is finite except 
at s = 0, in numerical implementation using finite difference, one can set the default intensity at 
the layer s = 0 to be the large value M instead of an infinite value. 
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4 Defaultable Equity Derivatives with Counterparty Risk 


We now apply our valuation methodology to value a number of defaultable equity claims. Specifi¬ 
cally, we will derive and compare the MtM values with and without counterparty risk provisions as 
well as the CRF value. Moreover, we will analyze and illustrate the bid-ask prices. 


As a special case of (2.2), we model the pre-default stock price process by 


where we assume constant interest rate r and default rates A®, i 6 {0,1,2}. In addition, we let 
A = Ylk =o ^ k ^ an d set 


a = L 2 A^(l - d 2 ) + - Li A^ (<J 2 - 1)+ + c 2 d 2 , 
P = L x A«(l - di)+ - L 2 A^ 2 ) (5i - 1)+ + ci <5i. 

4.1 Call Spreads 

Let us consider a generic call spread with the terminal payoff: 


(4.2) 

(4.3) 


g{S T ) = < 


m 2 if St > K + e 2 , 

{jn T^T 1 (S T -K) if K-e 1 <S T <K + e 2 , 
if S T <K-e i , 


ei+C2 

—mi 


(4.4) 


with mi, m 2 , ei, e 2 > 0, where m\/e\ = m 2 le 2 =: M. The payoff resembles that of a long position of 
M call options with strike K — e i, a short position of M call options with strike K + e 2 and short mi 
notional of zero coupon bond with the same maturity. Similar positions can be achieved when two 
OTC traders buy and sell call options with different strikes, plus/minus some cash. As ei and e 2 in 


(4.4) go to zero, the payoff converges to that of a digital call position covered in Henry-Labordere 


( 2012 ). 


With the terminal payoff g in (4.4), dividend h = 0, and value at reference default 1(tq) = 
—mi e~ r ^ T ~ T °\ the CRF value of the spread contract admits the formula 

n(t, s) =M (c BS (t, s;T,K — e\,r + A (0) , a) - C BS (t , s;T,K + e 2 ,r + A (0) , a)) - e _r ( T_t > mi , 

where C BS (t , s ; T, K , r, a) is the Black-Scholes call option price at time t with spot price s, maturity 


T, strike price K, risk-free rate r and volatility a. From (2.13), the MtM value with counterparty 
risk provision is given by 


P(t, s ) = E t ,; 


-(r+X )(T-t) 


g(S T ) + ^ e - ( r+A )(“ -t ) /(„, s u , P u ) du 


(4.5) 


where f(t,s,y) := A+ (A^ 1 ) + A^l — fi)y + (/5 — a)y + . The MtM value without counterparty 


risk provision P(t,s ) is similarly obtained replacing P u in (4.5) with IT 
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The model for S' in (4.1), the triple (g , h, l ), and other (constant) coefficients satisfy the conditions 
(C1)-(C7) with domain V = M + (see also (Heath and Schweizer 2000, Sect. 2)). We numerically 
compute the MtM value P(t,s) by Algorithm 1 from Section 3.2 For the iterative PDE (3.4), 
we adopt the coefficients in this section and the terminal payoff g(Sr) given in (4.4). In Table 
[2j we show the convergence of the MtM values with provision for three different contracts where 
d = 62 = {2,1,0.01}. The first column of each contract shows the value of the MtM value of the 
contract at spot s = 10 for each step 0 < n < 5. The second column of each contract shows the 
supremum norm e = ||p( n ) — T > ^ n_1 ^|| 0 o for each step 0 < n < 5. The algorithm stops at n = 5 for 
all three contracts. 



6l = 62 = 2 

6i = 6 2 = 1 

6i = 62 = 0.01 


Pi n i(0,10) e 

pi”) (0,10) 6 

Pi n )(0,10) 6 

n = 0 

0 

0 

0 

n = 1 

-0.1197 0.9048 

-0.1293 0.9048 

-0.1326 0.9048 

n = 2 

-0.1387 0.0992 

-0.1490 0.0992 

-0.1526 0.0992 

n = 3 

-0.1377 0.0060 

-0.1479 0.0060 

-0.1515 0.0060 

n = 4 

-0.1377 0.0002 

-0.1480 0.0002 

-0.1516 0.0002 

n = 5 

-0.1377 < 10" 5 

-0.1480 < 10" 5 

-0.1516 < 10" 5 


Table 2: Convergence of the MtM values with provision P(0, s) of call spread contract at spot price s = 10 
(at-the-money) and mi = m 2 = 1. Parameters: K = 10, T = 2, t = 0, r = 2%, a = 25%, X 1 ' 0 ' 1 = 3%, 
AW = 5%, A( 2 ) = 15%, Pi = 40%, R 2 = 40%, <5i = S 2 = 0, e = 10" 5 , S = 40, AS = 0.01, At = 1/1000. 


Let us visualize the convergence of the MtM value with CR provision p( n \0,s) in Figure [I] 
(left). Using the tolerance level e = ICC 5 for the maximum difference over each iteration, the 
algorithm stops after 4 iterations. As we can see, the price functions P^(0,s) and P^(0,s) over 
0 < s < S = 40 are not visibly distinguishable. 



Figure 1: (Left) Convergence of the MtM values with provision of a call spread P(0, s) for s £ [10,18]. 
(Right) Comparison of the three MtM values of a call spread {11(0, s), P(0, s), P(0, s)} over the spot price. 
Parameters are given in Table [2] 

In Figure [l] (right), we plot three different values n(0,s), P(0,s) and P(0,s) for 0 < s < S. 
As we can see, the ordering of these three values can change completely depending on the spot 
price. For example, for large spot prices, we observe that the CRF value dominates the other two 
MtM values, but it is lowest when the spot price is small. Furthermore, the value without provision 
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dominates the value with provision for high spot prices, and the opposite holds true for low spot 
prices. 

Next, we look at the sensitivity of the MtM values with respect to the counterparty’s or own 
default risk, collateralization ratio and effective collateral rate. In Figure [2j the MtM values are 
decreasing in the counterparty default rate (left) and increasing in the participant’s own default rate 
(right), as is intuitive. Note that the MtM value with provision moves more rapidly with respect 
to the counterparty default rate, but the MtM value without provision is more sensitive in the 
participant’s own default rate. 

In Figure [3] (left), an increase in 82 reduces counterparty-risk exposure, and therefore, increases 
the MtM values with and without counterparty risk provision. The rate of increase in contract value 
slows down when the collateralization ratio exceed 1. In the over-collateralized range [1,1.2], party 
1 is no longer exposed to the counterparty’s default risk. The increase in the contact value (from 
party l’s perspective) results from the possibility of collecting the excess collateral upon party l’s 
own default. 

In practice, if the participant’s funding cost rate is high, the effective collateral rate can be 
negative (see Burgard and Kjaer (2011)). This implies a net interest payment by the participant for 
the long position due to collateralization. As the effective collateral rate becomes more negative, 
the contract values with and without provision decrease as we observe on the right panel of Figure 

EJ 



Figure 2: (Left) The MtM values with and without provision are decreasing in the counterparty default 
rate with A^ 1 ) = 15%. (Right) The MtM values are increasing in A^ 1 ) with A^ 2 ) = 15%. The CRF value 
stays constant as A^ 2 ^ or A ( b varies. Parameters: = e 2 = 0.01, mi = m2 = 1, s = 15, K = 10, T = 2, 

t = 0, r = 2%, cr = 25%, A(°) = 3%, R l = R 2 = 40%, 6 x = S 2 = 0%, ci = c 2 = 1%, e = 10" 5 , AS = 0.05, 
At = 1/250. 


We illustrate the bid-ask prices P b and P s of a call spread in Figure |dj On the left panel 
where the participant is assumed to be default-free, we observe the dominance of the three prices: 
P s > n > P b . However, in the bilateral counterparty-risk case, the ordering of prices is different in 
in-the-money (ITM) and out-of-money (OTM) ranges. We see that n > P s > P b in the ITM range, 
but P s > P b > n in the OTM range. 
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Figure 3: (Left) The MtM values of a call spread with and without provision are decreasing in counterparty’s 
collateralization ratio <52- (Right) The MtM values of a call spread are increasing in the participant’s effective 
collateral rate ci £ [—5%,0%]. Parameters: e\ = e 2 = 0.01, mi = m2 = 1, s = 15, K = 10, T = 2, t = 0, 
r = 2%, er = 25%, A<°) = 3%, = 5%, A^ = 15%, R ± = R 2 = 40%, 5 1 £ {0% (left), 100% (right)}, 

62 = 100%, ci = c 2 = 1%, e = 10" 5 , AS = 0.05, At = 1/250. 



Figure 4: Call spread bid and ask prices with counterparty risk provision under (left) unilateral counterparty 
risk A^ 1 ) = 0%, A^ 2 ) = 15%, and (right) bilateral counterparty risk A^ 1 ) = 5%, A^ 2 ) = 15%. Other parameters 
are the same as in Figure [2] 


4.2 Equity Forwards 


Equity forward contracts are commonly traded in the OTC market. With stock S as the underlying 
asset, we consider a forward with maturity T. The initial forward price Fq is set so that the contract 
has zero value at inception. When the underlying stock defaults, the stock price goes to zero, and 
the buyer has to pay the discounted value e~ r ( T-t o).Fo at the default time. The contract cash flow 
is described by the triplet ( g,h,l ) = (St — Fq, 0, — e~ r As the underlying stock price 

fluctuates over time, the MtM value also varies. 


The MtM value of a long forward contract P(t,s) with provision (see (4.5)) is computed using 
Algorithm [l] In Table [3j we show the convergence result of the MtM value at time t = 1 when the 
stock price S\ = 20, with initial forward price Fq = 10. The first column of each case shows the 
value of the forward contract for each step n £ {1,..., 6}. The second column of each case shows 
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the error in terms of the supremum norm |jp( n ) — P^ n_1 ^||oo ° ver the whole domain [0,T] x T>, and 
with tolerance e = 10 ~ 5 the algorithm stops at n = 6 in both cases. The number of iterations may 
depend on the initial value threshold e and upper bound S. As we observe, for sufficiently 

large upper bounds S £ {30,40}, the convergent prices are the same. 


o 

CM 

II 

to 

to 

II 

CO 

o 

5 = 40 

n = 0 

Value e 

0 

Value e 

0 

n— 1 

8.6900 

17.5592 

8.6900 

26.4284 

n = 2 

7.6124 

2.1391 

7.6124 

3.2034 

n = 3 

7.6777 

0.1289 

7.6777 

0.1927 

n = 4 

7.6751 

0.0052 

7.6751 

0.0077 

n = 5 

7.6752 

0.0002 

7.6752 

0.0002 

n = 6 

7.6752 

< 10~ 5 

7.6752 

< 10~ 5 


Table 3: Convergence of the values of a forward contract when spot price Si = 20 at t = 1, with maximum 
stock price S £ {30 (left column),40 (right column)}. Other common parameters: F 0 = 10, T = 3, r = 2%, 
cr = 25%, A ( °l = 3%, AW = 5%, A< 2 > = 15%, R 1 = 40%, R 2 = 40%, 8 i = S 2 = 0, e = 10~ 5 , AS = 0.05, 
At = 1/500. 


At time t, when the stock price is s, the CRF value of a long forward is given by 

n(t,s) = (a- e" r(T_i) F 0 ) • (4.6) 

In order to compute the MTM value of a long forward contract P b without counterparty risk 


provision, we apply (4.6) to (2.14) and obtain 

rT 


p (t,s) = n(t, s) + e m 


e -(r+A)(«-t) _ a ) n+( u , s u )) -pu(u, s u ) du 


— n(t, s ) + e t t s 

rT 


f e -(r+A.)(«-*) ( (/3 _ a)(5u _ e - r ^F 0 ) + - P(S U - e- r ^F 0 )) du 
II(M) + J ((/? - a)Et, a [e" r (“-*)(£„ - - p(s - g-Kr-u)^)) du. 

To simplify the above equation, we notice that 


E 


t,S 


e- r{ - u - t \S u - e~ r( - T - u ^F 0 ) + } = C BS (t , s ; T, e^^Fo, r + A®, a ). 


We apply similar arguments to the seller’s MtM value, and summarize as follows. 

Proposition 4.1. The bid-ask prices without counterparty risk provision, P b (t,s) and P s (t,s), of 
a stock forward without counterparty risk provision are given by 

P b (t,s ) = U(t,s) + J ([p-a)C BS (u,s]T,e- r( ~ T - u '>F 0 ,r + \(°\a)-P(s-e- r( - T -VF 0 )^ du, 

P s (t, s) = II(f, s) + f ^(a — P) C BS (u, s ; T, e - r ( T -“) p 0 , r + A*- 0 - 1 , cr) — a{s — e -r ( T-t ) Fq)^J du, 


with II(t, s) satisfying (4.6). 
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The fair forward price makes the MtM value of the contract equal to zero at the start of the 
contract. The CRF value in (4.6) implies that Fq = e rT S q. However, the fair forward price in 
presence of bilateral counterparty risk is found implicitly. Precisely, the fair forward price Tq 
of a long (resp. short) position makes the MtM value with counterparty risk provision satisfies 
P b ( 0, So; F o) = 0 (resp. P*( 0 ,5 0 ; F§) = 0 ). 

In Figure [HJ we plot the bid-ask prices P b and P s of the forward contract with counterparty risk 
provision together with the CRF value. On the left panel, all three values increase as the underlying 
stock price increases. Similar to the call spread case, the price ordering changes from U > P s > P b 
in the ITM range to P s > P b > n in the OTM range. On the right panel, both MtM values decrease 
significantly as the counterparty default rate increases. However, the MtM with provision moves 
more rapidly, similar to Figure [2j 



Figure 5: (Left) The value of a stock forward at time 0 is increasing in the spot price. The fair forward price 
is the spot price at which the contract value is zero. (Right) The MtM forward contract values of a stock 
forward are decreasing in the counterparty’s default rate A^ 2 ), with F 0 = 10, t = 1 and T = 3. Parameters: 
r = 2%, a = 25%, X^ 0) = 3%, A (1 ) = 5%, R x = 40%, A< 2 > = 15%, R 2 = 40%, 6 ± = S 2 = 0, e = KT 5 , 
AS = 0.05, At = 1/500. 


Remark 4.2 (Total Return Swap). Total return swaps (TRS) are also traded over the counter as 
an alternative to equity stock forwards. The swap buyer will receive the increase in equity value at 
swap expiration date T, St — So. On the other hand, the buyer continuously pays a premium rate 
p > r until the expiration date and also pays the decrease in the equity value at the expiration date 
to the swap seller. The TRS is represented as the triplet (St — Sq, —p, — So e~ r P _T °^). 


4.3 Claims with Positive Payoffs 

Some derivatives have positive payoffs, i.e. the triplet g,h,l > 0, including calls, puts and digital 
options. For both conventions, the nonlinear property of the price function disappears since we can 
substitute the nonlinear functions P + (t,s) and n + (t,s) in (3.6) and (3.8) by the linear functions 
P(t, s ) and n(t, s ). From this observation, we derive the formulae for the bid-ask prices of derivatives 
with positive payoffs. 
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Proposition 4.3. For any claim with g,h,l > 0, the bid-ask prices of with counterparty risk provi¬ 
sion satisfy 


P b (t, s ) = E t , s 
P s (t,s)=E t , s 


e -(r+cH-A(°)) (T-t) + J e -(r+a+A(°)) (u-t) + A (0) 5u )) dti 

e -(r+^A ( «))(T-i) ff (5 T ) + [ T e -^+h+^ 0) ){n-t)^ s ^ + A (0) ^)) 


(4.7) 

(4.8) 


Without counterparty risk provision (see Remark 2.1), the bid-ask prices satisfy 


P b (t,s) = U(t,s)-E ty 


t 


ae -(r+\) (u-t) n( U) S u ) du 
p e -P+X) («-t) n (u,S u )du 


p s (t,s) = n(t,s) — E t)S 

We apply these results to call options on a defaultable stock. 

Example 4.4. The bid price of a European call option with provision is given by 

p\t, s ) = e"“ (T -*) C BS {t , s ; T, K, r + A (0) , a). 
Moreover, the bid price of a European call option without provision is given by 


(4.9) 

(4.10) 


(4.11) 


P\t,s)= 1- 


OL 


\W) + A( 2 ) 


(l- e - (A(1)+A(2))(T " t) ^ C BS (t,s-T,K,r + X {0 \a). (4. 


12 ) 


In both (4.11) and (4.12), the bid prices have two components: the CRF value of the call option, 


i.e. Black-Scholes price, and a multiplier that depends on the parameter a (see (4.2)). First, suppose 
that a = 0. This happens, for example, when the contract is perfectly collateralized and effective 
collateral rate is zero, i.e. 62 = C 2 = 0. In this case, both bid prices are identical to the CRF 
value II(f, s) = C BS (t, s;T,K,r + A*^ , a) . The parameter a. is positive when the contract is under¬ 
collateralized by party 2, i.e. 62 < 1. In this case, the call option buyer is exposed to counterparty 
default risk, and consequently, the bid prices are less than the CRF value II(t, s). On the other 
hand, a becomes negative when the contract is over-collateralized by party 2 and the collateral rate 
is negligible, i.e. 62 > 1 and C 2 = 0. The call buyer receives additional financial benefit from excess 
collateral since the buyer only returns a fraction 1 — L\ of the over-collateralized amount at the 
buyer’s own default. As a result, the bid prices are greater than the CRF value. 

Figure [6] illustrates the two bid prices (with and without provision) for various maturities and 
spot prices. We observe the dominance of the price without provision over the price with provision. 
Moreover, the difference of two prices increase as maturity or spot price increases. This wider level 
of difference is attributed to the fact that the difference of counterparty risk exposure increases as 
the maturity increases. 

Next, we derive a number of price dominance relationships for the bid-ask prices and the CRF 
value of claims with positive payoffs. 
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Stock Price 



Maturity 


Figure 6: The bid prices of a call option with and without counterparty risk provision are increasing in stock 
price s £ [5,15] (left) and in maturity T G [0,5] (right). Parameters: A^ = 5%, A^ = 10%, So = 10, T = 
1, t = 0, r = 2%, a = 25%, R 1 = R 2 = 40%, K = S 0 = 10 and 5 1 = 5 2 = 0. 


Proposition 4.5. Consider any claim with the triplet g,h,l > 0. If a, (3 > 0, then the following 
price dominance relationships hold: 

P b {t,s),P s (t,s) <U{t,s) and P b (t, s), P s (t, s) <U(t,s) . (4.13) 

Furthermore, if A^ + A^) > a, (3 > 0, then we have 

P b {t,s) <P b (t,s) <U(t,s) and P s (t, s) < P s (t, s) <U(t, s) . (4.14) 


We provide a proof in Appendix |A.3[ In market conditions where the effective collateral rates 
ci,C 2 are much lower than the counterparty/own default rates A^, A*- 2 ), the conditions for (4.13) 
and (4.14) are satisfied. In turn, since the buyer’s (resp. seller’s) MtM value with provision assumes 
P b (resp. P s ) as the liquidation value, which is lower than the liquidation value II in the case 


without provision according to (4.13), this implies the price dominance relationships in (4.14), as is 
presented in Figure [ 6 j 


5 Defaultable Fixed-Income Derivatives with Counterparty Risk 

We now consider defaultable fixed-income claims under bilateral counterparty risk setting. We model 
the stochastic factor X by two affine diffusions: Ornstein-Uhlenbeck (OU) and Cox-Ingersoll-Ross 
(CIR) processes. Under the OU model, the factor process X follows 

dX t = k{ 0 - X t ) dt + a dW t , 
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with positive constant parameters k, 9, a > 0 that represent the speed of mean reversion, long-term 
mean, and volatility, respectively. In the CIR model, the factor process X follows 


dX t = k{6 - X t ) dt + oJx t dW t , 


where the positive parameters {k, 6 , a} satisfy k , 9 , a > 0 and the Feller condition 2 k6 > a 2 so that 
X stays positive at all times a.s. Both models for X satisfy the corresponding conditions among 


(C1)-(C7) (Heath and Schweizer 2000 Sections 2.1-2.2), and we will consider here swaps whose 


payoff are bounded continuous. 

A generic fixed-income contract is described by the triplet (g(x), h(x), l(t , x)), with default times 
{Tj}i=o,i ,2 as in (2.1) and Markovian intensities of the form A^ = X^\t,Xt). From (2.13), the MtM 
value with counterparty risk provision is given by 


P{t,x) = E m 


e - f^r(n)+\(u,x u ))du + / e -ft U (r(v)+X(v,X v ))dv f ^ x ^p^ du 


(5.1) 


where f(t, x, y) = h{x) + A (°\t, x) l(t, x) + (A ^\t, x) + A ( 2 \t, x) — f3{t , x))y + (/3(f, x) — a(t, x)) y + . 
The parameters a(t, x) and /3(t, x) are defined similarly to the definitions (2.10) and ( |2.11| ). Noticing 
that P(t,x ) in ( |5.1[ ) is a special case of P(t,s,x ) in (2.13), we use Algorithm [I] to find the MtM 
value by solving iteratively the following PDE: 


Qp(n) 

dt 


+ £ x P {n) - ( r(t) + A (i, x)) P (n) + fit , x, = 0 . 


(5.2) 


for (t,x) G [0,T) x M, with p( n \T, x) = g(x) for x 6 1 (x G M + for the CIR case). The MtM 
value without counterparty risk provision P(t,x) is similarly obtained by replacing P(u,X u ) with 


n(u, X u ) in the right hand side of (5.1). 

For numerical examples, we assume the constant interest rate r and constant default rates A^ 
and A^ 2 ). The reference default intensity is modeled by A^ 0) = {X t ) + A X under the OU or CIR 
model for X. In other words, we impose an upper bound X and a lower bound X_ = 0 on the 
reference default intensity, and this allows us to apply our contraction mapping result in Theorem 
3.3 For implementation, we apply the FDM and Algorithm [I] in Section 3.2 to (5.2) to obtain the 
MtM value with counterparty risk provision. 


5.1 Credit Default Swaps (CDS) 

After the 2008 financial crisis, CDS contracts are traded in a new standardized way whereby the 
protection buyer pays at a fixed premium rate, along with a non-zero (positive/negative) upfront 
payment at the start of the contract. A long position of a CDS contract written on the default event 
{to < T} is summarized by the triplet notation (0, —p, 1). 

In Table [4j we show the convergence of the CDS MtM values. The first and second columns for 
each model show the values of the CDS at x = 2% and x = 8%, respectively. In addition, the third 
column shows the error in terms of the supremum norm e = ||P^ n ^ — F > ^ n_1 ^|| 0 o over the entire grid 
for each step n = 1,2,... ,7. The algorithm stops at n = 5 under the OU model, and at n = 7 under 
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the CIR model with the common tolerance level e = 10 5 



OU 

CIR 



Pi") (0, 2%) 

PT)(0,10%) 

e 

PT)(0,2%) 

Pl n )(0,10%) 

e 

n 

= 0 

0 

0 

- 

0 

- 


n 

= 1 

-0.1326 

0.05218 

0.9048 

-0.1130 

0.1569 

0.3950 

n 

= 2 

-0.1526 

0.04518 

0.0992 

-0.1131 

0.1176 

0.0911 

n 

= 3 

-0.1515 

0.04562 

0.0060 

-0.1130 

0.1248 

0.0159 

n 

= 4 

-0.1516 

0.04560 

0.0002 

-0.1130 

0.1238 

0.0022 

n 

= 5 

-0.1516 

0.04560 

< 10" 5 

-0.1130 

0.1239 

0.0003 

n 

= 6 

- 

- 

- 

-0.1130 

0.1239 

< 10" 4 

n 

= 7 

- 

- 

- 

-0.1130 

0.1239 

< 10" 5 


Table 4: Covergence of MtM values of a CDS in the OU/CIR model with reference default rate x = 2% and 
x = 10%. For the OU model: 9 = 3%, a = 3.5%, k = 5%. For the CIR model: 9 = 3%, a = 5%, n = 5%. 
Other common Parameters: T = 5, t = 0, p = 100 bps, r = 2%, ijj 0 = 0, i/; 1 = 5%, i^ 2 = 25%, w 0 = 1, 
Wl = w 2 = 0, Ad) = 5%, A (2 ) = 25%, R 1 =R 2 = 40%, ^ = 5 2 = 0, e = 10“ 5 , X = 20%, 1 = 0, AA = 0.001, 
At = 1/500. 




Figure 7: The CDS upfront prices under the CIR Model are increasing in the counterparty recovery rate R 2 
(left) and decreasing in the counterparty default rate A^ 2 ) (right). Other parameters are same as in Table [ 4 ] 
along with x = 8%. 


In Figure [7j we compare the three MtM values in terms of the counterparty recovery rate 
R 2 = 1 — L 2 and the counterparty default rate A^ 2 ). The CRF value is given by (A.9) and is 
independent of R - 2 and A*- 2 ), so it is flat across the x-axes in Figure [ tJ An increase in counterparty 
recovery rate or a decrease in counterparty default rate improves both MtM values, whereas the 
CRF value does not change. In both cases, the MtM value without provision dominates the MtM 
value with provision. 


5.2 Total Return Swaps (TRS) 

Total return swaps (TRS) on defaultable bonds are OTC traded, and their MtM values are subject 
to counterparty risk. A TRS is also referred to as a bond forward (see Schonbucher 2003 Chap. 2.5). 
Fix a maturity of T years, we consider a TRS on a zero-recovery defaultable bond with maturity 


T' > T. The value of the defaultable bond, denoted by C, is given in (A.5) for the OU model and in 
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(A.6) for the CIR model. Given no default up to the swap maturity T, the swap buyer receives the 


difference between the bond price C(T,Xt', T') and the pre-specihed strike K = C(0. Ay; T r ) from 
the swap seller. If the bond defaults before time T, the swap buyer pays the strike K to the seller. 
Until the first default time r or expiration date T, the buyer continues to pay at the risk-free rate 
plus a spread r + p. For the MtM upfront value P(t,x ) with counterparty risk provision, we follow 


(5.1) with the triplet 


g(x) = C(T,x ; T') — K, h(t,x ) = — K(r + p), l(t,x) = —K. 



Figure 8: (Left) CDS bid-ask upfront prices under the CIR model. (Right) TRS bid-ask upfront prices under 
the CIR model. Other common parameters: x = 8%, T = 3, T' = 10, t = 0, p = 100 bps, r = 2%, 6 = 3%, 
k = 5%, ipo = 0, ip! =5%, ip 2 = 25%, w 0 = l, Wl =w 2 = 0, AW = 5%, R ± = 40%, A^ = 25%, R 2 = 40%, 
Si = 82 = 0, e = 10 -8 , X = 20%, X = 0, AX = 0.001, At = 1/500. 


In Figure [8j we obtain the bid-ask upfront prices with provision for a CDS and a TRS, along 
with the counterparty risk-free upfront values, under the CIR model. Since both CDS and TRS are 
swaps, the MtM values can be positive or negative, depending on the reference default intensity. As 
buying a CDS is similar to longing default risk, the CDS upfront value is increasing in the reference 
default intensity (left). On the other hand, since the TRS buyer is shorting default risk, the TRS 
upfront value decreases as the reference asset’s default intensity A^ increases (right). 


6 Conclusion 

In summary, we have discussed a valuation framework to analytically study and numerically compute 
financial contracts subject to reference and counterparty default risks with bilateral collateralization. 
In addition to the underlying price dynamics used in this paper, our fixed point approach and the 
corresponding iterative numerical algorithm can potentially be adapted to price derivatives with 
counterparty risk in other models. The challenge lies in efficiently and accurately solving a sequence 
of inhomogeneous PDE problems. Our model also sheds light on the role played by counterparty 
risk and collateralization in the formation of bid-ask spreads. 

Several interesting and practically important problems remain for future investigation. For 
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example, what are the price and risk impacts of utilizing multiple counterparties for OTC trading? 


A recent study by Bo and Capponi (2013) derives the XVA of a CDS portfolio with a large number 


of entities. Apart from pricing, counterparty risk should also be incorporated into static or dynamic 
strategies for derivatives trading (see e.g. Jiao and Pham (2011)), as well as mark-to-market timing 


Leung and Liu (2012). In view of the financial crisis, counterparty risk has also become a key 


component in the design of clearing houses. Answers to these questions will be useful not only for 
individual or institutional investors, but also for regulators. 


A Appendix 


A.l Proof of Theorem 13.31 

For w G Cj’ 2 ([0, T) x V, M) and the operator M. in (|3.1[), we define v = v(t, s, x) by 


v := Mw = 


e~ $ g(S T ,X T ) + / e-^^ d v f ( u , Su , Xu ,w(u, SuiXu ))du 


(A.l) 


Equation (A.l) admits the same form as (1.2) of Heath and Schweizer (2000). By Theorem 1 of 
Heath and Schweizer (2000), v is a classical solution (C 1,2 ([0,T) x X>,M)) of the PDE: 


dv(t, s, x ) 

Ft, 


+ Cv(t, s, x ) — f(t, s, x) v(t, s, x) + f(t, s, x, w(t, s, x)) = 0 . 


(A.2) 


for (t, s, x) G [0, T) x D under certain conditions. To apply their result, we verify the sufficient 
conditions (A1),(A2),(A3') and (A3a') — (A3e') in their Theorem 1. The conditions (Al), (A2), (A3') 
and (A3a') — (A3c') are identical to our conditions (CT), (C2), (C4) and (C5) — (C'6). Since w G 
C^’ 2 ([0, T) x P,M), it is Lipschitz-continuous on [0, T] x D n , Vra. Combined with condition (C 7), 
it implies that the composition ( t,s,x ) — > f(t,s,x,w(t,s,x)) is uniformly Holder-continuous on 
[0, T] x D n x M, thus satisfying (A3cf). Lastly, the boundedness of v from Lemma 3.1 corresponds 
to (A3e'). Therefore, we conclude that v is a bounded classical solution (i.e. C^’ 2 ([0,T) x D,M)) of 
the PDE ( |A.2 ) for ( t , s, x) G [0, T) x D. 

Now, let’s select the initial function P G (7^’ 2 ([0,T] x e.g. = 0. Then, the 

subsequent functions P^ = AiP^ n ~ l \ n = 1,2,..., are also C}’ 2 ([0, T] x D,M) and satisfy the 


linear inhomogeneous PDE (|3.4|). By Proposition 3.2, the contraction mapping AJ ensures the 


sequence (P( n >) to converge to a unique fixed point P G C 5 QO, T) x V, 


A.2 Proof of Proposition [473] 

First, we denote A := A^ 1 ) + X < P and h(t, s ) := h(t, s ) + A^ 0 ^ l(t, s ). Applying the positive payoff to 
the definition of P b = P in (2.13), the MtM value with CR provision is given by 


P b (t,s)=E ty 


rT rT 

e -(r+A)(T-t) 5 ( 5r ) + / e -(r+A S u ) du + / (A - a) e -( r+A ) ( u ~Vp b {u, S u ) du 

(A.3) 


To prove (4.7), we substitute it into the RHS of (A.3) and verify that it indeed reduces to (4.7). 
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To this end, we get 


P b (t,s) = E M 


g -(r+A) (T-t) 5 ( 5t ) + J e -(r+X){u-t)~ h ^ Su)du 

+ / T (A - a) e- (r+A) (M - t} ( e -(r+a+X^){T-u) + F e ~(r+ a + A<°>) (v-u)^ gj dy \ du 


= E- 


t,S 


= E 


t,S 


= E t ,s 


e _ {r+Q+A (o)) (T-t) 5 ( (St ) 

pT pT pv 

Jt Jt Jt 

g —( r+a+A (o ))(T-t) 5 ( 5T ) 

/■T /-T 

+ J e-( r+x ^ u -^h(u,S u )du + J (e~( r+a+xm )( u -V-e-( r+x ^ u ~V)^u,S u )du 

e -(r+a+ A (0)) ( T — t ) 5 ( 5t ) + [ T e ~( r + a + A (0 >) 5 u ) du 


Since the last equality resembles (4.7), we conclude. The same steps will yield the proof for expression 
To verify (4.9), we use the expressions of II in (2.4) and P h = P in (2.14) to get 


(4.8). 


P 6 (M) =E m 


e -(r+A)(T-t) 5 ( 5t ) + / e -(r+A)(«-t)^ U)S g du _ / Q g -(r+A) (u-t) n ( u , $,) du 


+ / A e -( r+A ) n(u, 5 u )d- 


u 


= Et, 


0 -(r+A)(T- 


/■T /-T 

t) 5 (5 T )+ / e- (r+A)(u - f) h(u,5 M )du- / ae -(r+ A )(«-*)n( U ,S u )du 

Jt Jt 


+ Jt A e" (r+A) ( “" t} 


= E 


t,S 


f e -(r+A ( °))(«-u) 5„) du + e -(r+AW)(T-«) 

V ix 

e -(r+A(°))(T-t) _ jT T a e -(H-A) («-t) n(«, £ u ) du 

(r+A) + [ T (l _ e -M»-i)) e -HA(°))(«-4) h(u, 5 W ) dv 


du 


+J'P- 


= E 


£,S 


e -(r+AC<«)(T-t) 5 (g r ) _ a e -(r+A) (u-i) gj du + ^ e -(r+A<°>) (u-t)fc( Uj gj du 


= n(t,s) -E t) . 


ae -(r+A) («-i)n(u,S’„)du 


Applying the same steps to the definition of P s in (2.17), we obtain the equation (4.10). 
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A.3 Proof of Proposition |4.5| 

From (4.9) and (|4.10|) and the condition a, (3 > 0, we obtain the inequalities P b (t,s ) < II(t, s) and 


P s (t,s) < II(t, s). The price expressions (4.7) and (4.8) imply that 


P b (t,s) =E m 


.l t T (r+a+\W)du g ^ + f e -S^r+ a +\W)dv~ h ^ uSu)du 


< IE t,s 


e-ft ( r+x ^ du g(S T ) + [ e-ft U P+ x(0) ) d v~h( u , S u ) du 


= n (t,s) 


(A.4) 


Similar arguments give P s (t,s ) < II(i,s). Hence, we conclude (4.13). 

Next, applying the definition of P b = P in (2.14) along with the inequality (A.4), we get 


P b (t,s) = 


r T _ r T 

e -(T-t)(r+A) 5 (5 T ) + / e -(u-t)(r+A )h( u , S u ) du + / (A - a) e-( u ~V ( r+A In(u, S u ) du 

Jt it 

r r T r T 

>E M e - (T -' )(r+A) 5 (S r )+ / e-( u - t ^ r+X) h(u 1 S u )du+ / (A -a)e-( u ~ tnr+x) P b {u,S u )du 

l Jt Jt 

= P b (t,s). 


The last equality follows from the dehnition of P b in (2.13). From the definition of P s in (2.17) and 


the inequality 0 < P s < n in (4.13), we obtain 

rT 


P s (t, a) =E t>J 


> E tiJ 


-(T-t) (r+A) + j e -{u-t){r+X)~ h ^ Su)du+ f (A - 0) ( r + A )n(u, S u ) d U 

it Jt 


e -(T-t)(r+ x ) g (S T ) + / e -( u - t) ( r+x) h(u,S u )du+ / P)e~ {u ~ tHr+X) P s (u,S u )du 


= P S (t,s) 


The last equality holds from the definition of P s in (2.15). Hence, we conclude (4.14) 


A.4 CRF CDS and TRS Prices with OU and CIR Reference Default Intensities 

Here we summarize the computation of the CRF values of CDS and TRS. Let us assume that the 
risk free rate is a time-deterministic function r(t), and the default intensities i 6 {0,1,2}, are 
of the form AW(t, Xt) = ipi(t ) + Wi Xf , where X follows the OU or CIR model. This specification is 
more general than the one used in the numerical examples in Sections |5.1]]5.2| and it yields analytic 
expressions for the pre-default prices (without counterparty risk provision) of the credit default 
swaps and total return swaps considered therein. These prices are used for comparison analysis (see 
Figures [ 7 ] and [8]) . 

In the OU model, the pre-default zero-coupon bond price with maturity T and zero recovery 


meets the formula (Schonbucher 2003, Chap. 7.1.1) 


Ci (i, x] T) = e~ £ r(M) du e~ £ Mu)du E 


t,X 


e — f t wo Xs ds 


— e ~ ( r (u)+i/j 0 {u))du A\ (T—t)—Bi (T—t) wo x 


(A.5) 


where 


P'Uy 1 1 Av 'LL 

A 1 (u) = [-g 2 Bi{v) 2 - K0Bi(v))dv, Bi(u) = --—, 

./n ^ K 


0 < u < T. 
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In the CIR model, the bond price is given by (Schonbucher, 2003, Chap. 7.2) 


C 2 (t, x;T) = e- £ r(u) du e~ £ M^du^ 


g- ft W 0 X n du 


_ g - (r(u)+ip 0 (u))du _ t)e~ B ' 2 ^ T ~ t ' >X , 


(A.6) 


where 


A 2 {u) = 


2Eez 


(S + k) (e“' u - 1) + 2H 


R 2 (w) = 


2 (e s '“ - 1) mo 


.(E + k) (e s '“ — 1) + 2Hj ’ 


(A.7) 


for 0 < u < T with constant H = V^ 2 + 2 cr 2 mo ■ 

In turn, we can summarize the prices for the CDS discussed in Section 5.1 Under the OU model, 
the buyer pays the pre-default upfront price for the CDS: 

n (t,x)= FcA^x-, u) (w 0 xe~ K + (k 0- —)(ru-t)^^(l-^ K ^)-p] du , (A.8) 

Jt \ K K ) 

where C\ (t, x ; u) is the pre-default zero coupon bond price with maturity T and zero recovery 
given in (A.5). Under the CIR model, the pre-default upfront price of a CDS with maturity T and 
premium rate p is given by 


II(t, x) = / C- 2 (t,x;u) wq(k9 R 2 (u — t) + B 2 (u — t)x) — p 


du . 


(A.9) 


with C 2 (t, x ; u) in (A.6) and R 2 (-) in (A.7). See Chap. 7 of Schonbucher (2003). 

As for the total return swap described in Section 5.2, its CRF upfront value at time t < T, is 
given by 


II(f,x) = E m 


e -if (r(u)+AM(ti,X„))du/^^ Xt . ^ 


- [ e ~ f t u (r(v)+\(°Hv, x v))dv Kdu- I 

Jt Jt 

= C(t, x ; T') — K ^1 + pj C(t,x;u)du S j 


■ ft (r{v)+\W (v,X v )) dv K + du 


(A.10) 


For the CRF prices in Figures [ 7 ] and [8]), we set ip o = 0, and assume constant interest rate r. The 
counterparty default rates A^ 1 ) and A^ are also constant, and the reference default rate is set to be 
A(°)(t,Xf) = X t . This amounts to setting mi = m 2 = 0, and mo = 1 in (A.5)-(A. 10) to obtain the 
CRF values. 
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